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ABSTRACT 



Context. Optical interferometry is a powerfiil tool for observing the intensity structure and angular diameter of stars. When combined 
with spectroscopy and/or spectrophotometry, interferometry provides a powerful constraint for model stellar atmospheres. 
■ Aims. The purpose of this work is to test the robustness of the spherically symmetric version of the Atlas stellar atmosphere program, 

SAtlas, using interferometric and spectrophotometric observations. 

Methods. Cubes (three dimensional grids) of model stellar atmospheres, with dimensions of luminosity, mass, and radius, are com- 
puted to fit observations for three evolved giant stars, ijj Phoenicis, y Sagittae, and a Ceti. The best-fit parameters are compared with 
previous results. 

Results. The best-fit angular diameters and values of are consistent with predictions using Phoenix and plane-parallel Atlas mod- 
els. The predicted effective temperatures, using SAtlas, are about 100 to 200 K lower, and the predicted luminosities are also lower 
due to the differences in effective temperatures. 

Conclusions. It is shown that the SAtlas program is a robust tool for computing models of extended stellar atmospheres that are 
^ ' consistent with observations. The best-fit parameters are consistent with predictions using Phoenix models, and the fit to the interfer- 

1/^ ' ometric data for tf/ Phe differs slightly, although both agree within the uncertainty of the interferometric observations. 

■ 
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1. Introduction opacity distribution functions while the other uses opacity sam- 

ryQ pling; in this work we use the first version. The radiation field is 

Optical interferometry accurately measures the combination of computed using the meth od sugge sted by Rybickjl ( 11971ft . which 

^ the angular diameter and the structure of the intensity distribu- jg a reorganization of the iFeautrier ( 1964) method. 
K*" tion of a stellar disk. The combination of interferometry with 

k> spectroscopy and/or spectrophotometry is a powerful tool for de- I" *is work we compute cubes (three dimensional grids) of 

termining effective temperatures and other properties of stars, spherically symmetric stellar atmospheres with dimensions of 

and the growth of optical interferometry is testing theoretical luminosity, mass, and radius to model interferometric visibilities 

■ - - models of stellar atmospheres in new ways. to fit *e VLT WINCI observatiori s and b roadband spectropho- 

In a series of articles, IWittkowski et alJ ([200l 120063^ '"^r''^ froin |Johnson & Mi Mgia (ll975). This will provide a 

J rp , T~ — 7 ' — 777; rfT; — f robust test of the program and a comparison to the results pre- 

fit Very Large Telescope Interferometer (VLTl) observations ,. ,. -t^ ,,t. 

f 1 • . ■ \u \/TXT/-^T ■ . . -.u u ■ dieted using spherically symmetric Phoenix models. In the next 
of cool giants using the VINCI instrument with spheri- .. ° ., , ^ .. ., .. . ., ., 

cally symmetric stellar atmosphere models in local thermo- 



section, we outline the method for computing synthetic visibil- 



. . . /r rrr-\ .J -^u n ities and predicting radii, effective temperatures, luminosities, 

dynamic equilibrium (LTE) computed with the Phoenix pro- ^ , • . , , . 

JFt — 4. i l NnnnK j ..i. i ii i j masscs, and gravities. We present the results in the third section 

gram (IHauschildt et al.l 11999ft . and with plane parallel mod- . , .• • • ■ j- ■ 

T c n J A HP □ i m^nl hr!n^ ^ -ru ^nd the discussion in the fourth section. 

els from Phoenix and Atlas f Kuruca 119701 11993 ). The au- 
thors demonstrated that optical interferometry can detect the 
wavelength-dependent limb-darkening of cool giants and that 
the limb-darkening can be used to constrain model stellar atmo- 
spheres. These works are based on observations of three stars: 

t// Phoenicis (M4 III), y Sagittae (MO 111), and a Ceti (Ml. 5 111) 2. Method for Computing Visibilities and 
also called Menkar Determining Stellar Parameters 

The purpose of this note is to model these interferomet- 
ric observations using the new spherically symmetric version The goal of this note is to determine the global properties of 
of the Atlas program (SAtlas) developed by Lester & Neilson the three observed stars. Interferometric observations are ideal 
( l2008h and to determine fundamental parameters of the three for meas uring the angular diame ter, and with the HIPPARCOS 
stars. There are two versions of the SAtlas program, one using parallax dPerrvman & ESAII 19971) one can find the linear radii of 
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Table 1. Best-Fit Parameters of the Three Stars 





ijj Phe 


J Sge 


a Cet 


eLoCmas) 


8.7 ± 0.3 


6.3 ± 0.1 


12.6 ± 0.2 


^Ross(mas) 


Q 1 ^ -1- n 1 


o.uz ± u. i 


iZ. lU ± U.UJ 




1.66 


0.634 


0.98 


;r(mas) 


10.15 ±0.15 


11.9 ±0.71 


14.82 ±0.83 


R(R^) 


85 ± 1.6 


54 ±4 


88 ±5 


r.ff(K) 


3415 ± 87 


3650 ±78 


3584 ± 166 


L(Lo) 


882 ± 96 


468 ± 81 


1147 ±249 


M(Mo) 


0.85 ±0.1 


0.9 ± 0.2 


1.3 ±0.4 


logg (cm s--) 


Q 5^+0.06 




66+" " 



the stars. The visibiHty (iDavis et al ]l20Q(]HTango & Davisll2002h 
at wavelength A is 



(1) 



In this equation, 5^ is the sensitivity fun ction of the instrument 
(VINCI in this case) dKervella et alj200(]|) . Jq is the zeroth-order 
Bessel function, B is the length of the baseline of the interferom- 
eter, and ji is the cosine of the angle between the center of the 
stellar disk and a distance from the center The center of the disk 
corresponds to ji - \ and the edge of the disk is ju = as seen 
from the center of the star. To match the broadband VLTI/VINCI 
observations, the monochromatic model visibilities must be in- 
tegrated and squared to get the squared visibility amplitudes. 
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The denominator of Eq. |2] normalizes the visibilities, where 
is the flux of model. For a given model atmosphere, we can find 
the best-fit angular diameter with a minimum value. 

The next step is to use the model stellar atmospheres to fit 
spectrophotometric observations. The purpose of fitting spec- 
trophotometry is to constrain the effective temperature and the 
stellar flux; however, the effective temperature is dependent on 
the angular diameter. The modeled flux is scaled using the same 
method and assumptions as in Wittkowski et al. (2004) to fit the 
broadband data from Johnson & Mitchell (1975). The best-fit 
effective temperature is determined using the best-fit angular di- 
ameter from the interferometric fit to break the degeneracy. 

The linear radius of the star is determined using the best-fit 
angular diameter and HIPPARCOS parallax, and the luminosity 
is calculated using the radius and the effective temperature. The 
mass and gravity are constrained by comparing the predicted ef- 
fective temperatu re and lum inosity with theoretical stellar evo- 
lution tracks (iGirardi et al.ll200a) . 

3. Results 

For each of the three stars, we compute a cube of models with 
luminosity, mass, and radius as input parameters. The metallic- 
ity is assum ed to be solar, consistent with the conclusions of 
[Feast et al J (I1990) . and the microturbulence is zero. These as- 
sumptions can be tested, but for this work variations of metal- 
licity and microturbulence are ignored. The range of values for 
the input mass, luminosity and radius are chosen based on the 
results of Wittkowski et al. (2004, 2006a b). 

For 4i Phe, the mass range is chosen to be 0.6 to I.6M0 in 
steps of O.2M0; the luminosity range is 630 to ISSOL© in steps 




8.2 8.4 8.6 8.8 9 9.2 9.4 9.6 9.8 
Sld (mas) 




!.05 8.1 8.15 8.2 8.25 8.3 

Sross (mas) 



35 8.4 



Fig. 1. Values of minimum for each model stellar atmosphere 
for fitting the interferometric data of if/ Phe as a function of (Top) 
Sld and (Bottom) 6Iross- 



of 5OL0 and the radius range is 60 to 120/?© in steps of 2QRq. 
There are 480 models computed for ijj Phe. The mass range for 
y Sge is 1.0 to 1.9Mq in steps of O.IMq, the luminosity range 
is 400 to 7OOL0 in steps of SOL© and the radius range is 30 to 
7O/?0 in steps of IQRq, leading to 350 models for this star There 
are 1680 models computed for a Cet to provide a larger range of 
parameters and a more robust test of the SAtlas program. The 
mass range is 1.5 to 3.OM0 in steps of 0.1 M©, the luminosity 
range is 10^ to 2 x lO^L© in steps of 5OL0 and the radius range 
is 60 to lOO/^o in steps of IQRq. 

We fit the limb-darkened angular diameter, ^ld, to the inter- 
ferometric observations of each star and find a minimum value 
of The limb-darkened angular diameter is the angular diam- 
eter to the edge of the stellar disk corresponding to the the layer 
where fi - Q, and from the stellar atmosphere models we deter- 
mine the Rossland angular diameter by multiplying ^ld by the 
ratio of radius of the stellar atmosphere model at tross = 2/3 
to the radius of the outermost shell of that model, although this 
outermost shell is determined by an arbitrary choice of the mini- 
mum tross when the model is computed. The quality of fit to the 
interferometric data is sens itive to the Rossland angular diame- 
ter, as was noted by Wittk owski et al.l (|2004|) . The minimum 
values are shown in Fig. [1] as a function of limb-darkened and 
Rossland angular diameter for each model for the case of i/r Phe. 
The best-fit Rossland angular diameter is well constrained by in- 
terferometry and the values of the Rossland and limb-darkened 
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Fig. 2. (Top) The visibilities calculated for three S Atlas models 
and the interferometric data for t// Phe . (Bottom) Close-up of the 
second lobe of the visibility amplitude. The dotted line refers to 
the model stellar atmosphere for L = ISSOLq, M - 1.6Mq, and 
R - 6QRq, the dot-dashed line represents a model with L - 
980Lo, M = 1.2Mo, and R = SQRe, and the dashed Une is for 
L = 63OL0, M = O.6M0, and R = 1207^0. 



angular diameters, along with the uncertainty of the fit, for each 
star is given in Table [1] The Rossland angular diameter has a 
smaller uncertainty than the limb-darkened angular diameter for 
each star because the fits to the interferometric data produce 
less variation of the Rossland angular diameter than the limb- 
darkened angular diameter. This variation is clearly shown in 
Fig.IH 

The model visibilities are shown in Fig. |2] for iff Phe, with 
a close-up of the second lobe. Fig. [3] for y Sge, and Fig. |4] 
for a Ceti with a close-up of the second lobe shown. The dis- 
played model visibilities are chosen to be the smallest, middle 
and largest luminosity and gravity from each cube of models. 
The model visibilities for y Sge and a Cet agree with the results 
of|Wittkowski et al. ( 2006a b). For y Sge, there is not enough in- 
formation to test the limb-darkening of the model atmospheres. 
While the model visibilities using SAtlas and Phoenix both fit 
the observed interferometric data well within the uncertainties 
of the observations, the minimum value of the;^'^ from fitting the 
SAtlas models is smaller than that using the Phoenix models. 
This implies that there are small differences in the model atmo- 
sphere intensity structures predicted from each program. It is not 
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Fig. 3. The visibilities calculated for three SAtlas models and 
the interferometric data for y Sge. The dotted line refers to the 
model stellar atmosphere for L - JOOLq, M = 1 .9M0, and R = 
3QRq, the dot-dashed line represents a model with L = SSOL©, 
M - 1 AMq, and R - 50Rq, and the dashed line is for L - 4OOL0, 
M = I.OM0, andR = IQRq. 



obvious why the predictions using Phoenix models and SAtlas 
models differ in this case. 

Having fit the interferometric data, we next fit the broadband 
spectrophotometric data. One concern about fitting stellar atmo- 
sphere models to spectrophotometry is that the effective temper- 
ature and Rossland angular diameter are related. In Fig. |5] we 
compare the best-fit values of Rossland angular diameter as a 
function of the best-fit effective temperature for each model for 
fitting tfr Phe. The predicted angular diameter is almost constant 
with respect to effective temperature for the fit to interferometry. 
Therefore the interferometric fit of the Rossland angular diame- 
ter is used to constrain the predicted effective temperature. This 
process is repeated for the other two stars and the best-fit effec- 
tive temperatures are given in Table [1] 

The next step is to calculate the linear radius and luminosity 
of the stars. Using the HIPPARCOS parallax with the predicted 
Rossland angular diameter, we determine the radius of each star, 
the values of the parallax and radius are given in Table [1] The 
stellar luminosity is determined using the effective temperature 
and radius, L = AnR^crT^^. The best-fit effective temperature 
and luminosity for ea ch star is plo tted in Fig. |6] along with evo- 
lutionary tracks from lGirardi et al.l (i2000l) to determine the mass 
and gravity of each star. 



4. Conclusions 

The purpose of this work is to test the spherical version of the 
Atlas program. By fitting interferometric and spectrophotomet- 
ric observations using model stellar atmospheres, we derive the 
properties of tfr Phe, y Sge, and a Cet. The derived parameters 
are compared to earlier results. The angular diameters that are 
determined by fitting model stellar atmospheres to interferomet- 
ric observations are the same, except for a Cet, for which we 
predict an angular diameter that is 0.1 mas smaller The dif- 
ferences in predicted angular diameter are likely due to differ- 
ences in limb-darkening of the model atmospheres generated 
using SAtlas and Phoenix. Also our minimum value of the 
fit of the limb-darkened angular diameter for tfr Phe, 1.66, is 
smaller than the predicted value of;^^^ from the spherically sym- 
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Fig. 4. (Top) The visibilities calculated for three S Atlas models 
and the interferometric data for a Cet. (Bottom) Close-up of the 
second lobe of the visibility amplitude. The dotted line refers 
to the model stellar atmosphere for L - IOQQLq, M = 3.QMq, 
and R - 6QR0, the dot-dashed line represents a model with L - 
I5OOL0, M = 2.3Mo, and R = SO/Jq, and the dashed Hne is for 
L = lOOOLo, M = 1.5 Me, and R = IOORq. 



metric Phoenix models, 1.8. The S Atlas models for tfr Phe that 
best fit the interferometric data have eff'ective temperatures in 
the range of 4000 to 4 500 K, while the Phoenix models from 
IWittkowski etal] (120041) are all 3550 and 3600 K. The fit of the 
SAtlas models to spectrophotometric observations predict an ef- 
fective temperature of 3415 K conflicting with the prediction of a 
higher temperature. For the effective temperature range of 3400 
to 3800 K, the minimum values predicted by fitting SAtlas 
models are 1 .70 to 1 .72, similar to the values found using plane- 
parallel Atlas models. 

The fit of the SAtlas models to broadband spectrophoto- 
metric observations are used to determine the effective tempera- 
tures of the three stars because spectrophotometry is much more 
sensitive to the effective temperature than interferometry. The 
predicted effective temperatures are smaller than those found 
in the previous works, but this difference is due to the meth- 
ods used for determining the effective temperature, not differ- 
ences hi_flie_stenar_ataos2h programs. If we use the method 
from IWittkowski et al.l (12004 ) where the bolometric flux /boi is 



calculated by integrating the spectrophotometric data and then 
r^jj = 4/boi/(cr6'|^j,^), we would predict similar effective tem- 
peratures. Any variation in that case would be due to differ- 
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Fig. 5. The best-fit effective temperature as a function of 
Rossland angular diameter found by fitting spectrophotometry 
(-h's) and interferometry (x's) for t/r Phe. 





/I 1 1 1 1 






/ ,2.5 

/ / 1.4 1.2 1.0 0.8 


0.6 _ 


7 


/2.2 /y/^^/^yX4^ 






X2.0 //////y///// 






/^'V/// /// / / / / / 






//lAsy / / // 











3.6 3.58 3.56 3.54 3.52 3.5 3.48 3.46 
Log Te„ 

Fig. 6. Comparison of the derived effective temperatures and 
luminosities of the three stars relative the derived parameters 
from previou s analy ses. The Unes are evolutionary tracks from 
iGirardi et al.l (l2000l) . the points given by x's are the predicted 
effective temperatures and luminosities found using the SAtlas 
models and the plusses are the values found using Phoenix mod- 
els. 



ences in the calculation of the bolometric flux and differences 
in the angular diameter This suggests the effective tempera- 
tures here would differ by at most 1%. Our results also dif- 
fer because of the spectrophotometric data used. The differ- 
ence is smallest for i/r Phe because the Phoenix and SAtlas 
fits both use the same spectrophotometric data for the fitting 
and agree within the un certainty. For j S ge, Wittk owski et al.l 
(.2006b.) complement t he iJohnson & Mitchell (.19751) data wifli 
narrow-band data from lAlekseeva et aTT i 1997|) while for a C et. 
IWitfl'Cowski et al.' ('2006a') use optical (iGlushneva et al.lll998allbl) 
and infrared ([Cohen et al. 1996) spectrophotometry. The effec- 
tive temperature difference for the remaining two stars are about 
150 to 200 K. 

Because the angular diameters are similar, the predicted radii 
are also similar for both fits with SAtlas and Phoenix. Therefore 
the differences between predicted luminosities are due to dif- 
ferences in effective temperatures and thus due to the different 
methods for determining the effective temperatures. The lower 
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effective temperat ures and luminositi es imply lower masses 
when compared to iGirardi et alj (l2000h evolutionary tracks and 
smaller gravities. 

The SAtlas models are consistent with previous results. 
Fitting the models to interferometric and spectrophotometric ob- 
servations have provided a robust test of the spherically symmet- 
ric version of the Atlas program and have shown that the SAtlas 
program is a powerful tool for studies of stellar atmospheres. 
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